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ABSTRACT 


Decades of observations show that CMEs can deflect from a purely radial 
trajectory yet no consensus exists as to the cause of these deflections. Many of 
theories attribute the CME deflection to magnetic forces. We developed ForeCAT 


(Kay et al. 2013, 2015), a model for CME deflections based solely on magnetic 


forces, neglecting any reconnection effects. Here we compare ForeCAT predic¬ 
tions to the observed deflection of the 2008 December 12 CME and find that 
ForeCAT can accurately reproduce the observations. Multiple observations show 
that this CME deflected nearly 30° in latitude (Byrne et al. 1 2010; |Gui et ah 2011) 


and 4.4° in longitude (Gui et al. 2011). From the observations, we are able to 


constrain all of the ForeCAT input parameters (initial position, radial propaga¬ 
tion speed, and expansion) except the CME mass and the drag coefficient that 
affects the CME motion. By minimizing the reduced chi-squared, xh between 
the ForeCAT results and the observations we determine an acceptable mass range 
between 4.5xl0 14 and lxlO 15 g and the drag coefficient less than 1.4 with a best 
fit at 7.5xl0 14 g and 0 for the mass and drag coefficient. ForeCAT is sensitive to 
the magnetic background and we are also able to constrain the rate at which the 
quiet sun magnetic held falls to be similar or to or fall slightly slower than the 
Potential Field Source Surface model. 


Subject headings: Sun: coronal mass ejections (CMEs) 
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Introduction 


Predicting the deflections of coronal mass ejections (CMEs) is key for forecasting space 


weather effects. CME deflections can occur in the corona or interplanetary space (Kilpua 


et al. 2009; Isavnin et ah 2013 Rodriguez et ah 2011 Wang et ah 2002, 2004, 2006, 2014), 


or as the result of CME-CME interactions (Lugaz et ah 2012). Coronal deflections have 


long been associated with the presence of coronal holes (Cremades et ah 2006 Gopalswamy 


et ah 2009a; Mohamed et ah 2012), however, more generally, CMEs may be deflected by 


gradients in the background solar magnetic field (Gui et ah 2011| Shen et ah 2011) which 
tends to deflect CMEs towards the Heliospheric Current Sheet (HCS) on global scales. 


We develop a model, Forecasting a CME’s Altered Trajectory (ForeCAT) (Kay et ah 


2013, 2015), for coronal CME deflections determined by the background solar magnetic field 


which can, in general, produce CME deflections of comparable magnitude to observations. 


Kay et ah (2015) use an improved version of ForeCAT to show that magnetic deflections 


follow observed solar cycle trends: primarily latitudinal deflections occur near solar 
minimum but more variety in direction as well as larger magnitudes occur as the HCS 


becomes more inclined and the background magnetic field strength increases. In Kay et 


ah (2015), we find that a CME typically obtains a constant deflection velocity in the lower 


corona and continues deflecting through interplanetary space as the drag and magnetic 
forces cannot significantly reduce the CME’s motion when the drag coefficient is of order 
unity. Deflections in interplanetary space proportional to the inverse of the radial distance 
distance are the signature of a constant radial velocity and a constant deflection velocity 


and do not require additional interplanetary deflection forces (Kay et ah 2015). 


In this letter, we present the first comparison between ForeCAT model results and 
an observed CME. Observations determine most of the initial inputs for ForeCAT (initial 
location and shape, and propagation and expansion models), leaving only the CME mass 
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and solar wind drag coefficient as free parameters. We show that not only can ForeCAT 
reproduce the observed deflection, but it also constrains the unknown mass and drag 
coefficient as well as the background magnetic held. 


ForeCAT 


ForeCAT simulates CME deflections from the forces determined from a static 


background magnetic held (Kay et al. 2015). The hux rope of a CME is represented as 


a toroidal structure. To initialize a CME, ForeCAT requires the CME position (latitude, 
longitude, and tilt with respect to the solar equatorial plane), shape (height, width, and 
cross-sectional width), and mass. 

The most important free parameters relate to the propagation and expansion of the 
CME. ForeCAT uses a three-phase radial propagation model similar to that presented in 


Zhang & Dere (2006). The CME begins in a slow rise phase, followed by rapid acceleration, 


and finally a constant propagation phase. Kay et al. (2015) assume the simplest form of 


expansion, self-similar expansion, however, ForeCAT can utilize other expansion models 
based on constraints from observations. 

The forces responsible for the deflection, both magnetic pressure gradients and 
magnetic tension, are calculated at grid points across the face of the CME torus. The net 
force determines the deflection motion. ForeCAT includes a non-radial drag force which 


restricts the deflection. We use the standard hydrodynamic form of drag (Cargill et al. 


1996 Cargill 2004) with a drag coefficient proportional to tanh/3 (Forbes et al. 2006) and 


assume purely radial solar wind. The volumetric drag is calculated as 


_ 2C d tanh fipsw -> ^ 

f D —-- VCME ,nr\^CME, r i 

7 TO 


(i) 
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where Cj is the drag coefficient, b is the CME’s cross-sectional width, p S w is the solar 


wind density, and vcME,nr is the non-radial CME velocity. As in Kay et al. (2015), we 


ignore the factor of less than unity which describes the projection of the cross-sectional area 
onto the propagation direction. We use an empirical form for f3 determined from Figure 


1.22 of Aschwanden (2005): f3(R) = 2.515 (R — 1) L382 , with R being the heliocentric radial 
distance. 

ForeCAT uses simplified models for the background solar wind. The solar magnetic 


held is represented using a Potential Field Source Surface (PFSS) model (Altschuler & 


Newkirk 1969; Schatten et al. 1969 Altschuler et al. 1977) with the Parker interplanetary 


magnetic held included above the source surface. ForeCAT describes the solar wind density 


using a modified version of the Guhathakurta et al. (2006) model in which the density is 


the sum of a current sheet model and a coronal hole model, with the relative weight of the 
contributions determined by the distance from the HCS. We vary the background to explore 
the effect of different magnetic held prohles, which we discuss later. 


Observations 


While the 2008 Dec 12 has been extensively covered in the literature, we focus on the 
results of Byrne et al. ( |2010[ ) and |Gui et al.| (2011), hereafter B10 and Gil, as they offer the 
most complete coverage. Figure [l] shows the latitude, Stonyhurst longitude (Gil only), and 
radial velocity of the CME nose out to 50 R & for both the B10 data (blue squares) and the 
Gil data (red circles). Both authors utilize data from the Sun-Earth Connection Coronal 
and Heliospheric Investigation (SECCHI) on board the twin Solar Terrestrial Relations 
Observatory (STEREO) satellites. At the time of the eruption the satellites were separated 
by 86.7°. 
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Fig.0 includes error bars estimated from the scatter in the Gil and BIO observations. 
For the latitudinal data, the error bars represent the absolute difference between the two 
sets of measurements at each height. We find average deviations of 0.5° below 3 R 0 , 2° 
between 3 and 10 Rq, and 4° beyond 10 R Q . For the longitudinal data, the uniform error 
bars of 4° are determined from the scatter within the Gil data. 


4. Determination of ForeCAT Inputs 


ForeCAT CME’s are initiated at a height of 0.05 R & above photospheric polarity 
inversion lines (PILs). Using a finite number of harmonic coefficients causes a ringing 
effect in the PFSS model, and for coefficients up to order 90 this effect disappears by 1.05 
R & . Using an Helioseismic and Magnetogram Imager (HMI) synoptic magnetogram for 
Carrington Rotation (CR) 2077 we identify a PIL in the quiet sun (not associated with an 
active region) centered at 52° latitude and 81.7° longitude with a tilt of -13.8°, near the 
initial location of the observations. 


B10 find that the empirical relationship w(R) = 13 R 0 22 describes the CME’s angular 
width, w in degrees, as a function of radial distance, R in solar radii, which we use 
to determine the CME width, c. Both the CME height and cross-sectional width (the 


parameters a and b in Kay et al. (2015)) are set at constant fractions of the CME width 
(a = Ac and b = Be). The fraction A is estimated from the radial distance between 
the flanks and the nose of B10 and is set at one throughout the duration of the CME’s 
propagation. The cross-sectional radius of the CME is harder to constrain. We set the 
cross-sectional radius fraction B = - { and explore the effects of these chosen values in section 

ESI 


The bottom panel of Fig. [l] shows the observed radial velocity determined from 
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the change in position of the CME nose, which combines the effects of expansion and 
propagation. Since we know the expansion, matching ForeCAT’s nose velocity to these 
observations allows us to constrain the radial propagation model. ForeCAT initiates the 
CME with a radial propagation speed of 75 km s _1 . Once the CME reaches 1.75 Rq it 
begins constantly accelerating until it reaches a final propagation speed of 480 km s^ 1 at 18 
Rq. Beyond this distance, the CME moves with constant speed. 


5. ForeCAT Results 


We constrain all of the ForeCAT input parameters from the observations except for 
the CME mass and drag coefficient and vary these parameters to obtain a best fit with 


the observed deflection. Typical CME masses range between 10 14 and 10 16 g (Gopalswamy 


et al. 2009b) but observations suggest that the 2008 Dec 12 CME had a mass less that 


2xl0 15 g below 10 R Q (Carley et al. 2012 DeForest et al. 2013). The drag coefficient, 


Cd , is typically set near unity, but Cargill (2004) determine coefficients as high as 300 for 
individual simulations. Fig. [l] shows ForeCAT results for a CME with a mass of 7.5xl0 14 g 
and drag coefficient set to 0, which corresponds to the reduced chi-squared best fit to the 
Gil measurements, which we discuss in section |A4| 


Fig-0 shows that ForeCAT reproduces the latitudinal motion of the CME. We also fold 
decent agreement for the longitudinal motion beyond 5 Rq. The longitudinal motion below 
5 Rq does not match the Gil data, however, measurements of longitudinal deflections in 
the low corona are inherently highly uncertain as line-of-sight coronagraph obsevations 
integrate in the longitudinal direction, so we do not include these points when determining 
a best fit. 


We see a change in the direction of the longitudinal deflection below 2 R Q . As shown 
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(2015), the deflection is a combination of both local and global magnetic 


Kay et al. 

gradients. The local gradients initially cause the CME to move very briefly westward. 
This motion is small due to the weak local gradients of the solar minimum quiet sun. As 
the CME propagates out radially the local gradients are overcome by global gradients 
determined by the relative location of the coronal holes and the HCS. These global gradients 
are strong enough to slow down and change the CME’s deflection to a eastward direction. 


5.1. Variation with CME Mass 


Figure |2](a) and (d) show the effect of varying the CME mass on the latitudinal and 
longitudinal deflections in the same format as Fig. [l] The dashed lines represent a low 
mass CME case of 10 14 g and the dot-dashed lines represent a higher mass case of 2xl0 15 g. 
All three masses have similar latitudinal profiles with less massive CMEs deflecting more as 


the decrease in density allows for greater nonradial acceleration of the CME, as seen in Kay 


et al. (2015). The initial gradients deflect the low mass CME beyond the HCS but near 3 


R & the latitudinal deflection reverses direction, bringing it back toward the HCS. 


Fig. |2](d) shows very little difference in the longitudinal behavior of the high mass and 
best fit CMEs with the high mass CME deflecting slightly less both initially westward and 
back eastward. The low mass case again shows a significant increase in the total deflection 
but retains the same westward then eastward motion. 


5.2. Variation with Drag Coefficient 

Figures |2](b) and (e) show the effect of varying the drag coefficient (Ch in Eq. [l|. The 
drag coefficient is very small in the low corona when the plasma beta is high so the effects 
of scaling the drag coefficients are not noticeable for the first few solar radii. Beyond a few 
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Distance (R s ) 


Fig. 1.— Comparison of ForeCAT results with observations of latitude (top panel) and 
longitude (middle panel) versus distance for the 2008 December 12 CME. The bottom panel 
shows the radial velocity of the CME nose, which combines propagation and expansion, 
versus radial distance which is used to constrain the propagation model. The black line 


represents the best fit from ForeCAT, the blue squares represent the results of Byrne et al. 


(2010 

), and the red circles the results of 

Gui et al. 

(2011 
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Fig. 2.— Sensitivity of the latitudinal ((a)-(c)) and longitudinal ((d)-(f)) deflection to varia¬ 
tions in the CME mass ((a) and (d)), drag coefficient ((b) and (e)), and background magnetic 
field model ((c) and (f)). The solid line represents the ForeCAT best fit and blue and red 
points the observations, as in Fig. [lj In Panels (a) and (d) the dashed line corresponds to a 
mass of 10 14 g and the dot-dashed line corresponds to a CME mass of 2xl0 15 g. In Panels(b) 
and (e) the dot-dashed line corresponds to a drag coefficient of 2. In Panels [2|c) and (f)the 
dashed line corresponds to the scaled PFSS magnetic held model. 
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solar radii, we see that the increased drag causes a decrease in both the latitudinal and 
longitudinal deflection. The effect is more noticeable in the latitudinal direction as the 
CME deflects more in that direction. 


5.3. Magnetic Background 


Observations of Type II radio bursts suggest that the PFSS magnetic field may fall 


too rapidly with distance in the corona above active regions (Mann et al. 2003 Evans 


et al. 2008). Kay et al. (2015) show that ForeCAT is sensitive to the rate at which the 


background magnetic field decreases with distance. 


The 2008 December 12 CME occurs in the quiet sun, rather than above an active 
region, but we still explore the sensitivity of the deflection to the magnetic background. 
Fig. 0 shows the unsealed PFSS model results in a good fit with observations, however, 
a different magnetic background could potentially result in a better fit with a different 
parameters. Figure [2](c) and (f) compares the unsealed PFSS best fit (solid line) with the 
results for the same drag coefficient and CME mass, but using a magnetic background 
corresponding the the PFSS model increased by a factor of R below 2.5 R & . 


Both the latitude and longitude show that decreasing the rate at which the magnetic 
field decays with distance results in stronger deflections. For the CME mass and drag 
coefficient of the unsealed best fit, the deflections exceed the observed deflections. In section 


5.4, we explore whether a different set of CME mass and drag coefficient produces a good 


fit. 














12 


5.4. Determination of the Best Fit 

We determine the mass and drag coefficient of the best fit by computing the reduced 
chi-squared, xh which measures the variation between the ForeCAT model and the 
observations. We calculate xl as 

2 1 v ( Vobs-UFC ) 2 ^ 

JST (2) 

where N is the number of data, v is the degrees of freedom, y obs are the observed position, 
Dfc are the ForeCAT position, and a obs are the uncertainty as defined in section [3j 
Computing xl requires comparing y obs and ypc at the same radial distance so we linearly 
interpolate between the ForeCAT results to determine the ForeCAT values at the distances 
of the observations. 

Figure [ 3 ] shows contours of xl resulting from 400 ForeCAT simulations with masses 
between 10 14 g and 2xl0 15 g and drag coefficients between 0 and 3. In Fig. I a), the xl is 
computed using only the latitudinal points of Gil as they use the more commonly used 
technique to reconstruct the CME position. In Fig. [3^b), the xl is computed using only 
the COR2 longitudinal points of Gil. Fig. |3](c) uses both the latitudinal and longitudinal 
points to get the total xl used to determine the best fit. The white line corresponds to 
xl = 1.5. A xl ~ 1 indicates a good fit, but values significantly higher or lower than unity 
imply that the data is either under or overfitting the data. 

The latitudinal xl contours show an extended region which produces acceptable values 
of xl near unity. These contours show some degeneracy between mass and drag coefficient 
with lower masses requiring higher drag coefficients. 

The longitudinal xl values in Fig. [3^b) includes points as low as 0.12 due to the low 
number of longitudinal points and their large uncertainty. The longitude points alone 
cannot be used to constrain the CME input parameters. 
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Fig. 3.— Contours of the reduced chi-squared, xl, versus mass and drag coefficient. Panels 
(a)-(c) show results using the unsealed PFSS magnetic held model and Figs. [3][d)-(f ) show 
results using the scaled PFSS magnetic held model. Panels (a) and (d) show xl determined 
using only the latitude, Panels (b) and (e) show xl determined using only the longitude, and 
Panels (c) and (f) show xl determined using both latitude and longitude. Note the difference 
in the contour range for the longitude panels. The white line indicates where xl — 1-5 for 
all panels and in panels (d)-(f) the dashed white line corresponds to the region from the 
unsealed version. 
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When the latitudinal and longitudinal points are combined into a single measure of 
xl the region corresponding to xl — 1-5 increases in Fig. J3^c) due to the low xl of the 
longitudinal comparison. 


We use only the latitudinal xl to determine the range of acceptable parameters. By 
restricting the range of plausible input parameters to where xl < 1.5 we can restrict 
the CME mass to 7xl0 14 to 8xl0 14 g and a drag coefficient less than 1.4. The best fit 
corresponds to 7.5xl0 14 g for the mass and 0 for the drag coefficient. 


Figure [3^d) , (e), and (f), show the same as (a), (b), and (c), respectively, but for 
the scaled magnetic background. All panels have the same contour levels of xl as the 
corresponding unsealed panels. For low masses and small drag coefficients, the scaled 
background values of xl consistently exceed those of the unsealed background, however the 
majority of parameter space has a longitudinal xl below one. The scaled background does 
produce xl ~ 1 f° r the latitude for high masses, above 1.5xl0 15 g. This slightly exceeds the 


upper limit of 10 15 g below 6 R & determined by DeForest et al. (2013). However, it remains 
plausible that for the 2008 Dec 12 CME, which originated in the quiet sun, the magnetic 
field may decrease slightly less rapidly with distance than the PFSS model. 


5.5. Variation with Shape Parameters 

ForeCAT assumes that the CME height and cross-sectional width (parameters a and 
b ) maintain fixed ratios (A and B) with the CME width. The previous results used A = 1 
and B = j. We vary A and B and compare the resulting xl contours with the first case. 
Figure |4] shows contours of the xl f° r on ly latitudinal points (analogous to Fig. [3|a)) for 
four different cases. The solid white line shows the region corresponding to xl = 1-5 for 
that set [A,B] and the dashed white line corresponds to xl from the control case of [1, |] 
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from Fig. [3j These cases have longitudinal xt significantly below unity for most of mass 
and drag coefficient parameter space. 

Fig. [4] shows that changing the shape parameters either causes a slight difference in 
the range of xt near unity or results in a poor fit to the data such that xt never reaches 
unity. When the CME height is decreased relative to the width (Fig. [4](a)) the range of 
acceptable masses and drag coefficient decreases, and shows a shift toward higher masses. 
The decrease in height causes these CMEs to experience stronger magnetic forces so the 
range of acceptable parameters tends toward higher masses, as high as 10 15 g. The range 
of acceptable drag coefficient decreases. When the CME height is increased (Fig. |4](b)) , 
the range in parameter space yielding acceptable xt shifts toward smaller masses, with 
acceptable values as low as 5xl0 14 g. The increase in height causes these CMEs to initially 
experience weaker magnetic forces so only low mass cases can reproduce the extensive 
observed deflection. 

Decreasing the cross-sectional radius (Fig. |4](c)) causes a decrease in the range 
of acceptable parameters which is shifted toward lower mass, about 4.5xl0 14 g. If the 
cross-sectional width decreases further the range corresponding to xt = 1-5 continues 
to decrease and eventually disappears. Increasing the cross-sectional radius results in a 
significant shift of the region corresponding to xt = 1-5 to masses larger than 1.7xl0 15 g, 
suggesting that the observations cannot be reproduced with a cross-sectional width much 
larger than that used in Fig. [3|c). We find that despite not being able to accurately 
determine several CME shape parameters we can still constrain the CME mass and solar 
wind drag coefficient. 
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Fig. 4.— Contours of the latitudinal reduced chi-squared, xh versus mass and drag coeffi¬ 
cient (analogous to Fig. [3|a)). Each panel uses a different set of shape ratios A and B. The 
solid white line indicates where xt = 1-5 for that A and B and the dashed white line shows 
xl = 1.5 from the case in Fig. |3](c). 
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6. Conclusion 

The ForeCAT results reproduce the observations of the 2008 December 12 CME and 
by computing the reduced chi-squared, we can constrain the CME mass and drag 
coefficient. The CME shape is only constrained, rather than uniquely determined from 
observations. CME masses between 4.5xl0 14 and lxlO 15 g and drag coefficients between 
0 and 1.4 correspond to good fits with xt near unity. We find that the observations can 
also be reproduced with a magnetic background that decreases slightly less rapidly with 
distance than the PFSS model. This comparison shows promise that ForeCAT can be 
successfully used to predict CME deflections or constrain CME properties (mass), or probe 
the solar background (drag coefficient and magnetic field). Future development will be 
done to improve the expansion and propagation models so that they do not need to be 
predetermined from observations. 

The authors thank J. P. Byrne for providing the CME observations. This work was 
supported by grants from the Coordination for the Improvement of Higher Education 
Personnel (Capes - Brazil, 99999.016405/2012-09; to L.F.G.d.S.) The authors would like to 
thank the anonymous referee for his comments. 
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